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Qi Context. The photometric signal we receive from a star hosting a planet is modulated by the variation in the planet signal with its 

orbital phase. Such phase variations (or phase curves) are observed for transiting hot Jupiters with current instrumentation and have 
' also been measured for one transiting terrestrial planet (Kepler 10 b) and one nontransiting gas giant (Ups A b). Future telescopes 

, (JWST and EChO) will have the capability of measuring thermal phase curves of exoplanets, including hot rocky planets in transiting 

and nontransiting configurations and at different wavelengths. Short-period planets with a mass below 10 Re are indeed frequent, and 
I— I ' nearby targets (within 10 pc) are already known and more are to be found. 

' Aims. We test the possibility of using multiwavelength infrared phase curves to constrain the radius, the albedo, and the orbital 

, inclination of a nontransiting planet with no atmosphere and on a 1:1 spin orbit resonance. 

Methods. We modeled the thermal emission of a synchronous rocky planet with no atmosphere and its apparent variation with the 
orbital phase for a given orbital inclination. We assume that the planet is detected by radial velocity so its orbital period and minimum 
mass are known. We simulated observed noisy phase curves and then applied an optimization procedure to retrieve the radius and 
Q , albedo of the planet and the inclination of the orbit. 

» I , Results. Airless planets can be distinguished from planets having a dense atmosphere and their radius, albedo, and inclination (and 

' therefore true mass) can be retrieved from multiband observations with MIRI-JWST and EChO in the 5-15 lum range. The accuracy 

^ ' depends on stellar type, orbital distance, radius of the planet and inclination: hot and large planets on highly inclined orbit are favored. 

1 I ^ As inclination above 60° represents half of the randomly oriented orbits, the growing population of short-period, terrestrial-sized 

planets detected by radial velocity surveys and transits should offer several nearby promising targets for this method, including 
CN ■ planets GJ5 8 1 e,b, and HD40307b. 

^ ' Conclusions. Stellar activity is likely to limit the accuracy of this method, at least for some stars. It has not been taken into account 

I in this study, and its effects will have to be addressed in future works. 

00 '. 

^ I Key words. Methods: numerical, statistical - Infrared: planetary systems - techniques: spectroscopic 

rn ; 

2 1. Introduction with in 10 pc for K-G stars and 6 pc for M dwarfs (iBelu et al.1 
' '201 1). In the long run, nearby nontransiting planets and, in par- 

^ I Results from Kepler dHoward et"a Lll201 ll) reveal that 20-25 % of ticular, the potentially habitable ones will be studied by direct 

. MO through K dwarfs host a short-period planet (P < 50 days) detection w ith instrument deriv ed from pioneer projects suc h 

^ with a radius between 2 and 4 R®. For the same range of or- as Darwin ('Cockell et al.' '2009!), TPF-I ("Lawson et al.' '20071), 

^ bital periods, the same study shows that more than 18% of G-K TPF-C (Levine et al. 2009), or New Worlds (Cash et al. 200§, 

^ dwai-fs host a planet with a radius between 1 and 2 Re, which but hot planets are too close to their stai- for coronographs, 

_Cd ^ represents a lower limit, knowing that some transits ai'e not de- nuUers and occulters. However, we can study them by ob- 

tected for such small planets. Radial velocity surveys of F to serving in multiple spectral bands the variations in the unre- 

MO stars with HARPS have also unveiled a similar population: solved star-i-planet flux due to the changing phase of the planet, 

the debiased occurrence of exoplanets with Msin / < 30 de- assuming that intrinsic stellar variations can be distinguished 

creas es from ~ 45% at a period of ten days to ~ 10% at two from the planetai-y modulation. Phase curves of gas giants have 

days dMavoret a l.l l2011 ). At least two thirds of this population been observed for transiting configuration (see for instance 

consist of planets with Msin i < 10 Me. Knu tson et al.l2009l). but can also be n ieasured in a nontran siting 



For the hottest of these exoplanets (e.g. Kepler 10b and configuration (|Crossfield et alj|2pi0'). 'Cowan et alj (|2007l) mea 



Corot 7b), the geometric transit probabihty can reach 30 %, but sured phase variations for both transiting and nontransiting sys- 

it remains below 10 % for most of fliem. As a consequence, tems, and used these to constrain the inclination, albedo and re- 

the chai-acterization techniques requiring transit configuration circulation efficiency. Bathala et al. (201 1) have measured both 

miss more than 90% of these planets. This becomes a criti- the phase curve and the secondary eclipse of the terrestrial planet 

cal point when exploring the vicinity of the Sun. StatisticaUy, Kepler 10b. In this case, the Kepler's spectral band (0.4-0.9 A^m) 

transiting planets from this population are unUkely to be found means that the planetary flux measured is dominated by reflected 
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light while the present paper addresses the thermal phase curve. 
Nevertheless, it is interesting to note that it is already possible 
to measure a phase curve with an amplitude of about 5 ppm and 
also that the amplitude of the phase curve is clearly better char- 
acterized than the depth of the eclipse (see Fig. 13 of Bathala et 
al.), even though these two quantities are equal for a synchro- 
nized transiting planet with no atmosphere (for a nightside tem- 
perature of OK). The precision on the eclipse depth is indeed 
limited by the eclipse duration, which is about 9% of the orbital 
period. 

The phase- and wavelength- dependent emission of hot 
Jupiters has been simulated with various types of atmo spheric 
models (Barman et al. 2005; Burrows et al. 2010; Fort nev et alJ 
12006 ; Showman et al. 2009). The phase-dependent visible flux 
and its dependence on inclination has bee n studied f or nontran- 
siting hot Jupiters (Kane & Gelino 2011). Cowan & Agol (2011) 
studied the thermal phase curve of hot eccentric Jupiters with 
a semi-analytic model. The possibility of characterizing the at- 
mosphere of a nontransiting terrestrial exoplanet by its thermal 
phase curves has been studied by Selsis et al. (i201 1.) . 

In the present article we focus on the specific case of a rocky 
planet that has no atmosphere and is tidally locked in a 1 : 1 spin- 
orbit resonance on a circular orbit. Both radius and mass have 
been measured for some of the smallest known transiting plan- 
ets (Corot 7b, GJ1214 b, Kepler 10b, planets in the Kepler 11 
system). Obtained values already point to a broad diversity of 
objects in terms of mass-radius relationship, some probably be- 
ing volatile-rich (like GJ1214 b), while some (like Kepler 10b) 
seem to be dominated by denser material. It is reasonable to as- 
sume that large analogs of Mercury with no atmosphere exist 
within this population, either because they formed dry in the hot 
inner part of the protoplanetary disk or because they are too hot 
and receive too much X/EUV irradiation to keep an atmosphere. 
In the case of a synchronous planet on a circular orbit, the tem- 
perature map of the planet is constant and only depends on the 
surface bolometric albedo. The thermal emission measured by a 
distant observer can thus be robustly calculated and depends on 
the planet radius and albedo, the orbital inclination, as well as 
the type of the host star and the distance the planet is from the 
star For orbital inclinations over 60° (half of randomly oriented 
orbits), the amplitude of the phase curve is 87 % of the secondary 
eclipse depth. For orbital inclinations greater than 30° (87 % of 
randomly oriented orbits), the amplitude of the phase curve is 
50 % of the secondary eclipse depth. Therefore, the magnitude 
of the phase modulation is in most cases comparable to the oc- 
cultation depth but without the restriction of the occultation du- 
ration that limits the signal-to-noise ratio. 

The present work addresses the possibility of constraining 
the radius, the Bond albedo, and the orbit inclination using 
the thermal phase curve observed in multiple spectral bands 
by th e James Webb Space Telescope - JWST dOardner et alJ 
|2006|) or the Exop lanet Characterization Observatory - EChO 
dTinetti et al.ll201 11) . recently selected for the Assessment Phase 
for the Cosmic Vision program of ESA. 

2. Model 

We model the thermal emission of a planet with no atmosphere, 
on a circular orbit and tidally-locked in a 1:1 spin-orbit reso- 
nance. We assume a null obliquity, wh ich is consistent wit h the 
tidal evolution of short-period planets (iLeconte et al.ll201^ . We 
also assume a stable host star (see Sect. 15.41 for further discus- 
sions). Because insolation is constant at a given location on the 
planet, the surface temperature distribution only derives from the 



local radiative equilibrium and does not depend on the thermal 
inertia of the surface. Local temperatures are simply calculated 
by 

T(0)^[^^^^^^^^)\ (1) 

where is the stellar flux at the orbital distance of the planet 
(constant for a circular orbit), and the zenith angle of the point 
source star Considering an extended star ch anges the tempera- 
ture near the planet terminator (see Fig. 4 of iLeeer et al.ll20ll]) 
by producing a penumbra. The impact on the global emission 
of the planet and its apparent variations is, however, negligible 
even in the most extreme cases like Corot 7b, where the stellar 
angular diameter is 28°. This is because the global emission is 
dominated by the hottest substellar regions, which are not sig- 
nificantly affected. A is the bolometric surface albedo (assumed 
to be uniform and independent of the incidence angle, in which 
case it is equal to the Bond albedo of the planet), and cr is the 
Stefan-Boltzmann constant. We neglect geothermal flux, which 
implies that the night-side equilibrium temperature of the planet 
is zero. We can also set a geothermal flux O that adds up to the 
insolation and insures a minimum temperature r„„„ = (O/cr)'/"*. 
This, however, afifects the global thermal emission in an observ- 
able way only for unrealistically high values of <i>. 

The model calculates the infrared flux emitted by the planet 
in any direction, assuming a blackbody emission (the emissiv- 
ity of the surface is set to 1) and an isotropic distribution of 
the specific intensities. We use a 36 x 18 (typically) longitude- 
latitude grid, a resolution that is suflicient to produce smooth 
phase curves. The disk-integrated emission at a distance d is 
given by 

j 

where R is the radius of the planet, I,\{T j) is the blackbody spe- 
cific intensity (W vnT^ /im"' ster"') at the temperature Tj{9) and 
at the wavelength A, S j is the surface of the cell j, and aj is the 
angle between the normal to the cell and the direction toward 
the observer Only locations visible to the observer (cos aj > 0) 
contribute to the sum. Because the orbit is circular, the observing 
geometry is only defined by the orbit inclination /. 

Figure [T] presents thermal phase curves obtained with the 
model. We can see the influence of inclination, radius, and 
albedo for two wavelengths. In a nontransiting configuration, 
photometry will not yield the absolute value of the planetary 
flux (which can be inferred from secondary eclipses), but only 
its variation with the orbital phase. Figure |2] shows the peak am- 
plitude of the modulation as a func tion o f wavelength, called the 
variation spectrum by Selsis et al. (201 1"). Phase curves obtained 
in one single spectral band can be reproduced by various combi- 
nations of R, A, and /, but multiwavelengths measurements break 
this degeneracy. The disk-integrated thermal emission simply 
scales as R^ at all wavelengths, while A and / both have an ef- 
fect on the spectral distribution. The albedo directly controls the 
intrinsic spectral properties of the disk-integrated planetary flux 
through the distribution of surface temperatures, and thus deter- 
mines the wavelength of variation-spectrum maximum. The in- 
clination, on the other hand, controls the amplitude of the phase 
variations seen by the observer. This effect of the inclination is 
also wavelength-dependent, but it has a negligible effect on the 
position of the variation-spectrum maximum as seen in Fig. |2] 
Figure [3] presents this effect in a different way by showing how 
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changing the albedo or the inclination shifts wavelength of the 
emission maximum iA„ax) as a function of the orbital phase. 

It is therefore possible, in theory, to infer the radius, the 
albedo, and the inclination of a planet from the phase variations 
of its thermal emission, at least with noise-free observations. 

3. Methodology 

3.1. Producing a phase curve 

We consider a nontransiting exoplanet on a circular orbit (found 
for instance by radial velocity surveys). Some characteristics of 
the system are known to the observer: its distance, the stellar 
properties, the orbital period of the planet, and the ephemeris. 
The planetary radius, albedo, and the orbital inclination (as well 
as the true mass), on the other hand, are unconstrained. The ques- 
tion we address is the ability to retrieve these three characteris- 
tics using spatially-unresolved spectro-photometric observations 
covering one (or more) orbital period. With our model, we can 
produce these thermal phase curve variations and add unavoid- 
able noise to them. To calculate the stellar photon noise, we need 
to assume a collecting area, a data point binning time, a spectral 
resolution and a target distance (here set to 10 pc). We neglect 
here other stellar phenomena such as stellar spots, flares or pul- 
sations (see Sect. 15.4b . We assume that one complete orbit of 
the planet is sampled in 20 data bins of equal duration. We con- 
sidered the collecting areas of the James Webb Space Telescope 
(25 m^) and of EChO (1.1 m^). We consider ten 1 micron-wide 
spectral bins from 5 to 15 fj.m. We use th is for both the Mi d 
InfraRed Instrument (MIRI) on the JWST dWright et alJl2004 . 
and EChO spectrometers. We have similar results for the MIRI 
filters with three broad bands centered on 7.7, 10.0, and 12.8 yum, 
with widths of 2.2, 2.0 and 2.4 yum, respectively, corresponding 
to the filters of MIRI. This tells us that little spectral resolution 
is enough to break degeneracies. 

The considered MIRI spectrometer provides spectroscopic 
measurements over the wavelength range 5 - 28.3 yum, but the 
thermal noise of the telescope dominates above 15 fj.m. For the 
EChO spectrometers, the measurements could actually cover the 

0. 5-16 fj.m range, and integrating to wavelengths shorter than 
5 jum could improve the results, especially for hot planets (see 
Sect. |5]l, but we kept this 5-15 /urn range for comparison 
purposes. In addition, the planet/star contrast ratio drops be- 
low 10 ppm at wavelengths shorter than 5 fj.m for the coldest 
planet we consider (see Sect. |4|. A photometric precision bet- 
ter than 10 ppm, even for several days of integration may not be 
achievable for the forthcoming generation of infrared telescopes. 
Furthermore, the reflected light has to be taken into account at 
short wavelengths (below 3 fj.m) for the planets we consider. Our 
model uses the isotropic assumption for the planetary flux, which 
is a fair approximation for thermal emission but less valid for 
the scattered light (see Sect. 15. 11 1. Several additional parameters 
would need to be added to model the reflected signal correctly. 
The associated degeneracies and the wavelengthdependence of 
the albedo are such that extending the wavelengths range to 
reflected light is not likely to improve the retrieval. For given 
albedo, inclination, and radius, we produce n noisy phase curves 
for each spectral band. The only difference between these n re- 
alizations is the random stellar photon noise, and instrumental 
noise when included. For each noisy phase curve, we determine 
the values of R, A, and / that minimize using the downhill 
simplex method. From these n different estimates of R, A, and 

1, we determine their median value and the associated error de- 
fined as the smallest interval containing 95 % of the retrieved 



values (the "2-cr" confidence level if the values have a normal 
distribution). With zero eccentricity and obliquity, we only have 
these three parameters to constrain, but our method is applica- 
ble to more general configurations (not addressed in the present 
study), requiring the retrieval of more parameters. This is, for in- 
stance, the case with nonsynchronized oblique planets for which 
four additional unknown parameters control the thermal emis- 
sion: the rotation period, two angles for the rotation vector, and 
the surface thermal inertia. 

3.2. The instrumental noise 

In Belu et al. (2011) we showed that, on average, taking in- 
strumental noises into account divides the signal-to-noise ratio 
(S/N) by two compared with shot noise, for terrestrial-exoplanet 
transit spectroscopy and photometry with JWST. We therefore 
replicate here the modeling in this reference for both EChO and 
JWST/MIRI: read-out noise, dark current, optics thermal emis- 
sion, and reduction of exposure time due to readout rate. We do 
not take the zodiacal contribution into account, both domestic 
and at the target system, because its variations are not likely to 
be in the same regime of fr equencies as the sig nals considered 
here. Modeling jitter noise dPeming et al J 120091) is beyond the 
scope of the present work, but in-fligh t calibrations may v ery 
well be able to mitigate this effect (see iBallard et al.ll201 d, for 
latest methods and performances). The exact calculations are de- 
tailed very well in the program used to compute them, which is 
available on demand. 

Given the broad spectral range aimed at by EChO, it will 
require being divided into several sub-bands. Our 5 - 15yum 
band will be divided into two sub-bands (Table [T] Tinetti et al.l 
1201 lb Marc OUivier, priv. comm.). It has not yet been decided 
which detectors set (all-Si:As or MCT-i-Si:As) will be used for 
EChO. However, in our simulations, we consider only the option 
of all-Si:As detectors, since the other one gives comparable re- 
sults. The characteristics of JWST/MIRI are readout noise RMS 
19 e^/pixel, dark current 0.03 e"/s/pixel, and a full-well capacity 
10^ e-. 

4. Results 

The properties of the studied star+planet systems are summa- 
rized in Table|2l which gives the stellar masses, orbital distances 
(and periods) and corresponding equilibrium temperatures of 
the planets (ranging from 420 to 1640 K). The stellar masses 
represent the stars of the solar neighborhood (dominated by M 
and K stars. The range of orbital distances is chosen to offer a 
high S/N. We test our procedure to retrieve the radius, albedo, 
and inclination in different cases and with different instruments. 
Figure m shows the accuracy in retrieving the radius, the albedo, 
and the inclination of a planet for two instruments (EChO and 
JWST/MIRI), for different stellar types and for orbital distances. 
In all these plots, the X-axis gives the real values of R, A, and 
/, which are used to produce the noisy phase curves, and the Y- 
axis gives the median of the best-fit values. The error bar gives 
the 95 % confidence level (2 cr). A good retrieval should fall on 
the dotted line. For the albedo retrieval (left), the radius is fixed 
to 2Rq and the inclination to 60°, which is the median value for 
randomly oriented systems. For the inclination retrieval (mid- 
dle), the radius is fixed to 2/?® and the albedo to 0.1. For the 
radius retrieval (right), the albedo is fixed to 0. 1 and the incli- 
nation to 60°. All plots include stellar photon noise and instru- 
mental noise. The accuracy on the retrieved parameters depends 
both on the S/N that can be achieved in measuring the variations 
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in the planetary signal, but also in the shape of the spectral dis- 
tribution of the planet emission in the observation bands, which 
allows us to break the degeneracy between radius and temper- 
ature. Figure |5] shows how the S/N depends on the stellar mass 
and the orbital period. For this plot, the S/N is calculated as fol- 
lows. We integrate this variation spectrum over the 5-15 fim 
range and compare this signal with the photon noise of the star 
integrated over the same wavelength range. One can see that the 
S/N monotonically decreases as the orbital distance increases 
and monotonically increases with the stellar mass between 0. 1 
and 1 Mq. This trend generally explains the increase or decrease 
in the error bars when changing the orbital period or the stel- 
lar mass. This is, however, not always true because the ability 
to retrieve A, R, and ; also requires measuring the wavelength 
dependency of the spectra. An important factor is thus the posi- 
tion of the emission peak compared with the observation spectral 
window. This explains the turnover around IMq. An appropriate 
window can compensate for lower S/N and vice versa. As an ex- 
treme example, we consider a very hot planet observed at long 
wavelengths where the emission of the substellar area (domi- 
nating the phase curve) is in the Rayleigh- Jeans regime. In this 
case, despite a possibly very high S/N, the weak wavelength- 
dependency would result in an enhanced degeneracy between 
radius and temperature. This effect can be seen, for instance, on 
the albedo and inclination retrieval in Fig.|4^: error bars can be 
similar to or smaller for the less massive star or at larger orbital 
distance, contrary to what S/N suggests. This is because cooler 
planets have an emission peak closer to or inside the 5-15 fim 
domain and exhibit a better temperature signature. The position 
of the variation-spectrum peak can be seen in Figs.|6]and|7] This 
tells us that optimal retrieval will be achieved when the observa- 
tion window includes the shortest wavelength not affected by re- 
flection. In this article we did not adapt the window to the planet 
properties, so our retrieval is thus not optimal. 

For inclinations close to 90°, the flux becomes roughly pro- 
portional to sin /, hence weakly sensitive to /. For this reason, in- 
clinations in the 80 - 90° range cannot always be distinguished, 
and the uncertainty on the inclination increases with the inclina- 
tion despite larger modulations. 

Although the accuracy obtained with EChO is lower than 
JWST for a given number of orbits, EChO will have the pos- 
sibility of dedicating more time to a given target, reducing the 
S/N to levels comparable to or better than those achievable with 
JWST, hopefully with a significantly improved stability and the 
capability of observing spectral phase curves over a wider spec- 
tral domain simultaneously. Using the 2.5 - 5 fim domain would 
better characterize the hottest planets. On JWST, measurements 
in this window imply the use of NIRSpec, an instrument that 
cannot be used simultaneously with MIRI, while EChO spec- 
trometer should cover the whole 0.5 - 16 fim window. 

The results presented in this work assume that the observa- 
tion covers two orbital periods, which will certainly be required 
in practice to extract the periodic planetary signal from the stel- 
lar variability. We tested on several cases that the error on the 
albedo, the inclination and the radius decreases as 1 / V^, where 
is the number of orbits observed. Also for observations cover- 
ing orbits, the values of R, A, and / obtained for each individ- 
ual orbit (or for a number of orbits smaller than A^) provide an 
information on their dispersion. 



4.1. Constraining botli mass and radius 

If we combine the projected mass Msin/ measured from ra- 
dial velocity observations with the constraint on the inclination 
from the phase curve, we obtain an estimate of the true mass. 
Retrieved mass and radius and their associated uncertainty can 
be compared with theoretical models to assess the composition 
of the planet. Figure [8] shows mass-radius relations of ice/rock 
and rock/iron planets from Fort nev et al.l (12007ft and the range 
of mass and radii obtained from phase curves for two planets 
(R = 1.5 Rs), Msini = 4 Me and = 2 R^, Msin/ = 8 M®) 
around a 0.5 and a 0.8 Mq star with an orbital period of three 
days. The planets around the 0.8 Mq star logically give the 
best estimate of the composition. The uncertainty on the com- 
position is dominated by the uncertainty on the mass, and one 
should also include the error on the measurement of M sin / it- 
self. Interestingly, the uncertainty on the planet radius does not 
come directly from the uncertainty on the stellar radius as it 
does when the planetary radius is inferred from the primary tran- 
sit depth, in the case of transiting planets. Indeed, the depth of 
the transit provides the ratio RplR-k and an estimate of Rp there- 
fore requires a value for which comes with its uncertainties. 
When deduced from the thermal phase curve, the radius estimate 
is also affected by uncertainties (for instance on the luminosity 
of the star or on the orbital distance), but not coming directly 
from the radius. Errors on the stellar luminosity and the stellar 
mass (used to convert an orbital period into an orbital distance) 
do, however, result in retrieval errors. 

As another illustration, we show in Fig. |9] the result for 
three known exoplanets (HD40307b, GJ581 b, and GJ581 e, see 
Sect. 15.21 ). assuming that they are made of silicates (an assump- 
tion that is probably not realistic for GJ581 b, whose volatile 
content must be high considering its large mass). 

4.2. Effect of an atmospliere 

We do not know a priori whether the observed planet has an at- 
mosphere of not. As shown by Selsis et al. (201 1), the presence 
of an atmosphere can be inferred from the variation spectrum. 
We did, however, test our R - A - i retrieval procedure (which 
assumes no atmosphere) on phase curves computed for a planet 
with a dense atmosphere. The planet is a 1.8 R® rocky planet 
with a 1 bar CO2 atmosphere, in a 1:1 spin-orbit resonance on 
an 8-day orbit around an M3 dwarf, with 60° inclination. The 
structure of the atmosphere and the associated radiative trans- 
fer has been modeled with a 3 D GCM (global climate model) 
described in lSelsis etakl (1201 lb . 

Our model obviously fails to fit the phase curves variations at 
all wavelengths, which by itself shows that the "no atmosphere" 
assumption is wrong. If we treat each spectral band indepen- 
dently of the others, our model can obtain a set of R, A, and / 
with reasonable at some wavelengths but the retrieved values 
strongly vary from one wavelength to another Figure [TO] shows 
the retrieved values as a function of wavelength (no noise is con- 
sidered in this retriev al). As in the variation spectrum described 
by Selsis et al. (201 V), the signatures of molecular absorption can 
be seen in these plots, in particular the 2.7, 4.3, and 15 //m bands 
of CO2 in the "inclination spectrum". This is because different 
wavelengths probe different altitudes in the atmospheres, with 
different day-night temperature contrasts. In absorption bands, 
the flux comes from high altitudes with smooth day-night con- 
trast, while in atmospheric windows, the flux comes from the 
surface that exhibits a temperature distribution that differs less 
from the airless case. 
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As an example, we show in Fig. [TT] the results of the fit for 
three wavelengths. We only fit the variations and not the abso- 
lute flux, which is why the solution found by our procedure can 
produce a phase curve that is shifted vertically. Because we only 
consider synchronized planets, which is equivalent to a planet 
with no thermal inertia, our modeled phase curves cannot ex- 
hibit the phase shift that atmospheric circulation can produce. 
If we no longer consider synchronized planets and let the rota- 
tion rate and the thermal inertia be free parameters, our airless 
model may fit a monochromatic phase curve obtained with an at- 
mosphere with better agreement, by reproducing the phase shift. 
In this case, the displacement of the hot spot (compared with 
the exact substellar location) is due not to horizontal circulation, 
as in an atmosphere, but to vertical heat diffusion. However, re- 
trieved rotation rate and thermal inertia would also depend on 
the wavelength, as for the inclination, radius, and albedo. Also, 
for planets hot enough to produce an observable infrared phase 
curve, on circular orbits, synchronization occurs on extremely 
short timescales so there is no justification for considering non- 
synchronized solutions. A phase shift can thus be attributed to an 
atmosphere alone. This point will have to be stressed more care- 
fully for cooler planets, which are less subject to tidal forces, and 
which need to be observed by future direct detection techniques. 

If no photometric variations are seen despite accurate pho- 
tometric measurements, we may be able to infer that the planet 
has a dense atmosphere, as attempted on GJ876 b by Seager and 
Deming (2009). The absence of modulation can also be due to 
an inclination close to 0°, but in some cases inclinations lower 
than a given value can be rejected. This can be done using the 
measurement of the projected rotation of the star and assuming 
that planetary orbits remain close to the stellar equator (known to 
be wrong for many hot Jupiters), but also using d ynamical con- 
strain ts, which was done for the system s GJ876 (ICorreia et al.l 
120101) and GJ581 jMavoret al.ll2009ah . A lack of modulation 
can also be come from an extremely high albedo, which may be 
checked at short wavelengths. 

5. Discussion 

5.1. Validity of tine radiative model and reflected light 

The ability to constrain R, A, and / from spectral phase variation 
comes from the simplicity of the model that assumes an isotropic 
distribution of the thermal emission, uniform surface properties, 
and an emissivity independent of the wavelength and equal to 1 . 
We now discuss the validity of these assumptions. 

An emissivity value e lower than 1 but independent of A 
would only mean that the "effective albedo" we infer is in fact 
equal to(A-i-e-l)/e and does not affect the radius and in- 
clination retrieval. Variation in the emissivity with wavelength 
(for instance from the 10 fim silicate band), if significant enough 
to affect the planetary emission in an observable way, would be 
seen in the variation spectrum just as an absorption feature due to 
an atmosphere (see Selsis et al. 2011). A correction could thus be 
done a posteriori (probably yielding larger uncertainties on the 
retrieval). The variation in the surface albedo with wavelength 
influences our model only in terms of emissivity as our model is 
only sensitive to the bolometric surface albedo (the fraction of 
reflected energy integrated over the whole spectrum). We tested 
our retrieval algorithm that assumes a uniform albedo on phase 
curves computed with a nonuniform albedo. To each individual 
cell of the surface, we gave a random albedo using a normal dis- 
tribution. We did the test for a mean albedo of 0.1, 0.3, ...,0.9 with 
a standard deviation of 0. 1 . Error bars on the retrieved albedo are 



centered on the mean value of the normal distribution, and the er- 
ror bars are twice as broad as with a uniform albedo, but there 
is no noticeable effect on the inclination and radius retrieval. We 
did not test the effect of having large regions with different albe- 
dos. An interesting case to be addressed in the future is to con- 
sider a surface composition (and thus a surface albedo) changing 
at a given temperature (and thus incidence angle). This could be 
relevant for very hot planets like Corot 7b or Kepler 10b, which 
have a strong temperature gradient and whic h may have subste l- 
lar lava-oceans covering a < < 45° area dLeger et al.ll20Tll) . 

Because of its very low thermal inertia and the length of its 
solar day (~ 27 days), the temperature distribution on the day 
side of the Moon is similar to that of a synchronized planet, and 
brightness temperatures on the day side of the Moon follow the 
cos(0)4 law, where is the incidence angle (Lawson et al.l2000h . 
This means that the isotropic assumption for the emission is a 
good approximation (with no significant effect of roughness and 
craters), but also that variations on albedo with location are small 
enough. Light reflected by the moon does not, however, follow a 
Lambertian distribution (the reason the full Moon seems so "flat" 
to us). As other rough planetary surfaces with a low albedo, the 
diffuse reflectio n is described better by the Lommel-Seeliger law 
(lFairbairnl2005 ). Because of the non-Lambertian behavior of the 
reflected light (which can also include a specular component for 
high albedo), we restricted our work to wavelengths (A > 5fim) 
where the reflected component is orders of magnitude lower than 
the emission. However, the temperature distribution, and thus 
the thermal emission, can be affected by the dependence of the 
surface albedo on the incidence angle 9 found in the Lommel- 
Seeliger phase function. We tested the impact of this effect on the 
thermal phase curve and found it to be insignificant on the phase 
curve. The reason is that surface temperature is affected only at 
high incidence angle, near the terminator where the temperature 
is too low to contribute significantly to the disk-integrated emis- 
sion. Small departures from an isotropic behavior of the thermal 
emission do exist on the Moon and other planetary surfaces (due 
for instance to roughness and craters) but we neglect these effect 
in this study. Knowing that the retrieval of R, A, and / is feasi- 
ble, it will become necessary to address this question in further 
works. 

5.2. Existing candidates 

Among the published planets detected by radial velocity and at 
the time of writing, at least seven are potentially teiTestrial ob- 
jects with characteristics allowing a measurement of the phase 
variations with JWST or EChCQ: a minimum mass below 20 Mm, 
an orbital distance within 0.05 AU, and me mbership in a system 
closer than 15 pc. According to estimates dMavor et al.ll201 iT) . 
more such planets remain to be found. Expected detections with 
15 pc critically depend on the occurrence of planets around M 
stars, which is yet poorly constrai ned B The occurren ce rate 
as inferred from Kepler candidates (iHoward et al.ll20lH) shows 
no decrease (or even a slight increase) from K to MO dwarfs. 
Assuming the same frequency of exoplanets for M and K stars 
yields hundreds of candidates within 15 pc. Within the current 
sample of seven, three have zero ecc entricity (at the prec ision 
of the measur ements): GJ5 81 e and b (iMavor et alJl2009al) . and 
HD40307b ( Mayor et a l."2009b). The other four have eccentric 
orbits: GJ674b: e ~ 0.2 (Bonfils et al. 2007), GJ876d: e ~ 0.2 



' Unless their ecliptic latitude does not allow pointing. 
^ A study by iBonfils et alj i2011l) based on HARPS measurements, 
actually shows that 35 % of M stars planets host a 1-10 Me planet. 



5 



A.S. Maurin et al.: Phase curves of non-transiting terrestrial exoplanets. 2. Airless planets 



(lRiveraetal.lf2010h . 61Vir b: e ~ 0.1 (IVogtetal.ll2010l) and 
55Cnc e: e ~ 0.25 ('Fisch er et aDl2008l) . To model the phase 
curve of an eccentric planet requires including heat diffusion 
on the subsurface, which is controlled by an additional param- 
eter, the surface thermal inertia. It is the only parameter to add 
if we assume that th e planet rotates at its equilibrium rotation 
(iLeconte et al.ll20f o'). The planet can, however, be trapped in a 
spin-orbit resonance as is Mercury. We will address the case of 
noncircular orbits and nonsynchronized planets in a forthcoming 
work. 

The case of GJ581 is interesting because of its proximity 
(only 6.2 pc) but also because two planets (e and b) in this sys- 
tem could produce a detectable modulation. Planet 'e' is hot- 
ter and smaller (Msin i-1.9 M®, a - 0.028 AU) while planet 
'b' is cooler but larger (Msin/=16 M®, a - 0.04 AU) and 
more likely to possess a dense atmosphere. Planet HD40307b 
(Msin /=4.1 Me, a = 0.047 AU, M. = O.SMq, d ^ U pc) 
represents another fine target. 

5.3. Transiting planets 

For planets in transit configurations, the radius and inclination 
are akeady known. The albedo can be constrained from the spec- 
trum of the planet at the secondary eclipse. However, as is the 
case of Kepler 10 b (Batalha et al. 2011), the phase curve may 
be better determined than the secondary eclipse depth due to the 
short duration of the eclipse. Phase-curve fitting could then pro- 
vide additional constraints. The phase curve of Kepler 10 b can- 
not be fitted with thermal emission alone. At the wavelengths 
observed by Kepler (0.4 - 0.9 jjm), the reflected light dominates, 
and the large amplitude of the phase curve points to either a very 
large Bond albedo (~0.6) or to a highly anisotropic reflection. 

5.4. Stellar variability 

Stellar variability comes from surface brightness inhomo- 
geneities (mainly spots and bright plages), their evolution in 
time, and their modulation by the stellar differential rotation. 
These inhomogeneities are produced by magnetic activity and 
convection, and they produce wavelength-dependent noise and 
systematic variations that will make the extraction of the plan- 
etary phase curve (and thus the constraints on the planet prop- 
erties less accurate). Using Kepler visible li ghtcurves of ~ 
150, 000 stars over 33.5 days, B asri et alj (1201 ID find that more 
than 70% of them exhibit relative variations larger than 10"^. In 
the mid-infrared, preliminary studies show t hat this amplitude 
is expected to be lower by a factor of ~ 5 dBean et al. 2010). 
To measure the planetary phase with sufficient accuracy, peri- 
odic residuals of the stellar variability at the orbital period of 
the planets have to be lower than the planetary signal (which is 
larger than 10 ppm in this work). Period filtering is thus required 
to lower the amplitude of the noise by at least a factor of ~ 100. 

The efficiency of this filtering will depend on the variabil- 
ity power spectrum of the host star, its rotation period compared 
with the planet orbital period and the level of characterization 
of this variability that can be achieved. Evaluating the feasibil- 
ity of this extraction is beyond our expertise and the scope of 
this paper, but we stress here that it certainly represents the main 
challenge for this type of high-precision spectro-photometric ob- 
servations covering two or more orbital periods (days to weeks). 
It is obviously not an imp ossible task, at least for some quiet 
stars. Batalha et al. (1201 1|) were for instance able to extract the 



phase curve of Kepler 10 b, which has an amplitude of less than 
10 ppm in the visible range. 

6. Conclusion 

In this study we have modeled the thermal emission of an exo- 
planet with no atmosphere, on a circular synchronous orbit. The 
modulation of this emission with the orbital phase at different 
wavelengths can be observed for exoplanets in non-transiting 
configurations, which are known from a previous radial veloc- 
ity detection. We showed that these infrared multiband phase 
curves can be used to infer the absence/presence of an atmo- 
sphere and, in the airless case, to constrain the radius and albedo 
of the planet, as well as its orbital inclination. The constraint on 
the inclination yields a constraint on the mass using the value of 
M sin / measured by radial velocity. The knowledge of both the 
radius and the mass can then be used to assess the bulk compo- 
sition of the planet. 

Typical precisions of 10° on the inclination, 0.1 on the 
albedo, and 10 % on the radius can be obtained with JWST, and 
15° on the inclination, 0.2 on the albedo, and 15 % on the radius 
with EChO, for nontransiting large terrestrial planets with a pe- 
riod of a few days or less around nearby 10 pc K and M stars. 
This accuracy is obtained from observing two orbits (for an ide- 
ally quiet star) and can be reduced by by observing 2 
orbits. Assuming that a whole orbit is observed, the accuracy 
is not very sensitive to the sampling of the orbit. We used ten 
bands with a 1 //m width from 5 to 15 fim and the MIRI's filters, 
although the number of bands is not critical either. Two bands is 
the minimum to break the albedo-radius degeneracy. Better re- 
sults can be obtained for the hottest planets by including shorter 
wavelengths. 

These planets represent an abundant population, as shown 
by radial velocity surveys and Kepler More than 90 % of these 
hot terrestrial exoplanets do not transit their host star and may 
only be characterized by their thermal phase curve in the fore- 
seeable future. This method can be applied with EChO or JWST 
on three akeady known exoplanets (GJ58 1 e and b, HD40307 b), 
and more should be discovered by radial velocity surveys. For 
each transiting exoplanet with an observable secondary eclipse, 
there will be statistically ten times more similar objects closer to 
the Sun and with an observable thermal phase modulation. 

The retrieval of R, A, and / works obviously better for large 
and hot planets that provide higher planet/star contrast ratios and 
do so for highly inclined orbits that produce a larger amplitude of 
the phase variations. Thanks to the distribution function of ran- 
domly inclined orbits that is equal to sin /, planets are rare whose 
phase variations are not observable because of a low inclination 
(only 13% with/ < 30°). 

Unlike secondary eclipse observations, the exposure is not 
limited by transit duration. Therefore, phase curve fitting can 
provide better characterize of the planet properties than the 
eclipse, as in the case of Kepler 10 b. This method implies the 
ability to achieve spectrophotometry with 10"^ relative accuracy. 
This precision must be stable over the duration of the observa- 
tion used to extract the phase curve. Stellar variability must be 
characterized and removed in such a way that the residuals do 
not contaminate the planetary signature. Stability and stellar ac- 
tivity are the most critical issues for phase curve observations. 

Further modeling is required to include the effect of surface 
roughness and craters on the surface temperature end thermal 
emission. The case of eccentric planets will be addressed in a 
forthcoming paper. 
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Fig. 1. Effect of inclination (left, R - 2R^ and A - 0.1), radius (middle, / = 60° and A - 0.1) and albedo (right, / - 60° and 
R - 2 R^) on the flux ratio, at 8 (default) and 15 //m. Phase corresponds to the minimum fraction of the illuminated planet received 
by the observer (0 if / - 90°). The eclipse normally occurring for / = 90° is not shown here. 
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Fig. 2. Effect of inclination and albedo on the variation spectrum. 
The peak amplitude of the phase curves is given as a function of 
wavelength for a planet at 0.04 AU from a 0.5 star. Solid 
curves are obtained with an albedo of 0. 1 and inclinations rang- 
ing from to 80°. Dotted curves are calculated for an inclination 
of 60° and albedo ranging from 0. 1 to 0.9. 
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Fig. 3. Influence of inclination and albedo on the wavelength 
'l/iiflA of the spectrum peak. Calculations are done for a planet 
at 0.04 AU from a 0.5 Mq star, for inclination between and 80° 
and an albedo of 0. 1 (solid lines) and for an albedo between 
and 0.9 and an inclination of 60° (dashed). 
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5-11 fim 


11-15 jum 


Telescope 


Type of detectors 


RON 


dark current 


full-well capacity 


RON 


dark current 


full-well capacity 


EChO 


Si:As 


12 


0.1 


2-10^ 


19 


0.1 


2-10^ 


EChO 


MCT+Si:As 


1000 


500 


37-10" 


12 


0.1 


2-10^ 



Table 1. Characteristics of detectors used for EChO. The readout tioise RMS is iti e /pixel, the dark curretit iti e /s/pixel, atid 
full- well capacity in e". Only all-Si: As are used in the results shown in this paper, MCT-i-Si:As giving comparable results. 
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0.2 


0.5 


0.8 


0.02 


T=596 K 
P=2.31 d 


983 K 
1.46 d 


1641 K 
1.15 d 


0.04 


T=421 K 
P=6.53 d 


695 K 
4.13 d 


1160K 
3.26 d 



Table 2. Equilibrium temperatures (K) and orbital periods (days) 
for different stellar masses and orbital distances. Temperatures 
are calculated with a Bond albedo of 0.1. 
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20 40 60 



80 
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Fig. 4. Ability to retrieve the albedo, inclination and radius of a synchronous rocky exoplanet with a 1.2 m telescope like EChO, 
whose instrumental noise is simulated with Si:As detectors (rows: a, b) and with the MIRI/JWST (c, d) for different stellar masses. 
The orbital distance is a=0.02 AU (a, c) and 0=0.04 AU (b, d). The triangles, squares, and crosses correspond to 3 different stars 
(M=0.2, 0.5, and 0.8 Mq respectively). Error bar include 95 % of retrieved values. The median of the retrieved values is indicated 
inside the error bar 
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Fig. 5. Estimated signal-to-photon-noise ratio (S/N) as a func- 
tion of stellar mass and orbital period. Each curve in the left plot 
is for one orbital distance, and each curve on the right plot is for 
one stellar mass. Here, we consider the planetary signal to be the 
peak amplitude of the phase curve, integrated from 5 to 15 ;um 
intervals. The duration of the two data bins used to compute the 
peak amplitude is set to 1/20''' of the orbital period, the inclina- 
tion to 60°, the distance to 10 pc, and the telescope diameter to 
1.5 m telescope. 
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Fig. 6. Variation spectrum as planet/star contrast (left) and S/N (right). These plots give the peak amplitude of the phase variation as 
a function of wavelength. The orbital distance of the planet is 0.02 AU, its radius 2 R®, and its inclination 60°. Each curve is for a 
given stellar mass. To calculate the S/N, we set A/l = 1 yum, distance=10 pc, telescope diameter=1.5 m. The gray area corresponds 
to the 5 - 15/vm range (left & right) and to contrasts higher than 10"^ (left). 
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Fig. 7. Same as Fig. |6]for an orbital distance of 0.04 AU. 
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Ms=0.5 Mq Ms=0.8 Mq 



=40° 
= 60° 
=80° 
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M/Me 



15 
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Fig. 8. Constraining the bulk composition. Boxes indicate the 95 % confidence level on the retrieval of both the radius and the mass 
for two synchronous planets (R - 1.5 Rg), Msini - 4 M^, and R - 2 R(g, Msin; = 8 Me), at 0.02 AU of a 0.5 and a 0.8 Mq 
star, observed with the JWST for 3 different inclinations (2 orbits for computation). The solid, dashed and dotted lines correspond 
respectively to a pure iron, pure rock, and pure-ice planet according to Fortney et al. (.2007.) . 
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Fig. 9. Precision of the mass and radius determination for three 
known exoplanets (GJ581 b, GJ581 e, and HD40307 b). Similar 
to Fig. [8] but we no longer assume a fixed radius for the planets, 
but a rocky composition (100 % silicates). The radius changes 
with the mass and thus with the inclination. 
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Fig. 10. Best fit parameters as a function of wavelength, for a planet with a 1 bar CO2 atmosphere. 
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Fig. 11. Phase curves of a planet with a 1 bar CO2 atmosphere (solid lines) and best fit of the variation (and not of the absolute flux) 
obtained assuming an airless planet (dashed lines). The fit is done for each wavelength separately. 
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